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Abstract 

A systematic numerical investigation of a recently developed nuclear structure ap- 
proach is presented which diagonalizes the Hamiltonian in the space of the sym- 
metry-projected Hartree-Fock-Bogoliubov (HFB) vacuum and symmetry-projected 
quasiparticle excitations with respect to it. The underlying HFB transformation, 
which is assumed to be time-reversal and axially symmetric, is determined by vari- 
ation after the projection. The model allows the use of large basis systems. It has 
been applied to the calculation of energy spectra of several even-even, odd-odd and 
odd mass nuclei in the sd shell with mass numbers reaching from A = 20 to 30. 
The Chung- Wildenthal interaction has been used. Good agreement with the exact 
shell model diagonalization and a considerable improvement on a previous approach, 
where the HFB transformation was significantly more restricted, is obtained. 



1 Introduction 



Microscopic nuclear structure models try to explain the excitation energies and other 
experimental spectroscopic information starting from a realistic effective nucleon- 
nucleon interaction. An example is the interacting shell model. It is generally 
acknowledged to be the most fundamental theory of a nucleus [0] 0. But since 
the number of A-nucleon configurations to be diagonalized increases rapidly with 
the number of single particle states, its application is limited to rather small basis 
systems, like the sd shell . Thus for many problems the configuration spaces have 
to be truncated drastically. How to select a numerically feasible number of physical 
relevant configurations is the central question of all nuclear structure models which 
use the shell model as the conceptual basis. 

One prescription for such a truncation is provided by the mean field approaches like 
the Hartree-Fock |P or the Hartree-Fock-Bogoliubov (HFB) [0 § theory. They 
use the variational principle to describe the ground state of a nucleus by one single 
Slater determinant. However, usually in this determinant all physical symmetries 
are broken. Thus one has to use projection techniques to find the component with 
the desired quantum numbers in order to obtain physical states. Moreover, in order 
to obtain the energetically lowest solution for a given set of quantum numbers, 
one has to perform the projection on the good symmetries before the mean field is 
determined by variation. 

Following this idea, a whole hierarchy of symmetry conserving mean field theories 
and their extensions into multi-configuration mixing methods have been proposed by 
some of our group a few years ago 0. They are known as the VAMPIR (Variation 
After Mean field Projection In Realistic model spaces) and the MONSTER (MOdel 
for handling many Number- and Spin-projected Two-quasiparticle Excitations with 
Realistic interactions and model spaces) approaches. 



The models of the VAMPIR family, like VAMPIR |Tg, EXCITED VAMPIR |n 
FED (FEw Determinant) VAMPIR and EXCITED FED VAMPIR 0, can be 
used to describe the lowest few states of a given symmetry, irrespective of their 
particular structure. Each state is obtained by variational calculations. So, e.g., the 
VAMPIR model approximates the lowest state of certain quantum numbers by just 
a single symmetry-projected HFB vacuum. The underlying HFB transformation 
is determined by variation after the projection on the desired quantum numbers. 
Correlating configurations as well as excited states can be obtained by chains of 
similar variational calculations. 



On the other hand, the MONSTER approaches are constructed for the description 
of complete excitation spectra with respect to a one-body transition operator 0| 
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1^. Thus they are well suited, e.g., for the description of giant resonances. States 
which are reached by a one-body operator are of similar structure as the initial 
ground or yrast state. Therefore they can be described by expanding the nuclear 
wavefunction around a symmetry- projected reference vacuum, which may be either 
a usual HFB, or, e.g., a VAMPIR solution for the ground or yrast state. This is 
done in the MONSTER approach. More explicitly, the MONSTER configuration 
space is built out of the symmetry-projected vacuum and all symmetry-projected 
two-quasiparticle excitations defined on it. The excited states are then obtained by 
diagonalizing the residual interaction between these configurations. 



In all realistic applications up to now, these models have been simplified by im- 
posing symmetry restrictions on the underlying HFB transformations. In the first 
VAMPIR calculations only real, time-reversal invariant and axially symmetric HFB 
transformations, which neither mix proton and neutron states nor states of dif- 



ferent parity, were admitted [|ri|]. With this real VAMPIR approach, as it is 
called, only states in even-even nuclei with even spin and positive parity can be 
described. If a MONSTER calculation is based on such a real VAMPIR transforma- 
tion, states with different symmetries are described entirely by the two-quasiparticle 
components of the wavefunction. Some years ago, this VAMPIR approach has been 
improved by using HFB transformations, which are essentially complex and allow 
proton-neutron and parity mixing |T^. Only time reversal and axiality are still im- 
posed on the transformation. This complex VAMPIR approach can describe states 
of arbitrary spin parity in even-even and odd-odd nuclei. It considers many more 
nucleon correlations ||15||. 



Recently the MONSTER approach has been generalized for the use of such com- 
plex VAMPIR transformations [0. It was tested for two light even-even sd shell 
nuclei. The results were encouraging. The performance of the complex MONSTER 
description for heavier even-even nuclei, for odd-odd and odd-mass nuclei, however, 
remained an open question. In the present paper we try to answer this question by 
a systematic investigation. 

For this purpose the complex MONSTER model was applied to several even-even, 
odd-odd and odd mass systems. The calculations were performed again in an sd 
shell basis, in order to be able to compare with the exact shell model results. The 
considered nuclei of each type are covering the range from very light systems to 
the ones with the biggest possible numbers of shell model configurations in the sd 
shell. The results are compared to those of the older more restricted real MON- 
STER approach, too. It is demonstrated that the complex MONSTER approach 
approximates the shell model results very well and shows a clear improvement on 
the older real approach. 

In the next section the main features of the MONSTER on VAMPIR model are 
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summarized and the consequences of the symmetry restrictions for the HFB trans- 
formation are shortly explained. In section (|^) the results of the application are 
presented and discussed. Finally conclusions and an outlook are given in (^. 



2 The Model 



The detailed formulation of the MONSTER on VAMPIR approach has been given 
elsewhere [^. Therefore only the main ideas of the model are sketched here. 



2.1 The VAMPIR model 



A finite model space is defined by a D-dimensional set of orthonormal single par- 
ticle wave functions, which are eigenstates of some spherically symmetric one-body 
potential, e.g. the harmonic oscillator. In second quantization the corresponding 
creation and annihilation operators are denoted by {c|, cl, . . .}d and {q, Ck, ■ ■ ■}d, 
respectively. The indices i, k summarize the quantum numbers of a state. It is as- 
sumed that the effective many body Hamiltonian appropriate for this model space 
is known and can be represented by a sum of only one- and two-body terms 

H = f + V. (1) 

The one-body part T contains the matrix elements of the kinetic energy (or, if an 
inert core is assumed, the single particle energies) and the two-body part V the 
matrix elements of the effective interaction. 



Starting from the chosen particle basis, the Hartree-Fock-Bogoliubov (HFB) trans- 
formation 1^ is used to define the corresponding quasiparticle basis with its quasi- 
particle creation and annihilation operators, {aj^, ajj, . . .}d and {oq,, ajs, . . .}d- 
enumerate the quasiparticle states. The HFB transformation is the most general 
linear transformation which ensures that the Fermion anti-commutation relations 
are fulfilled by the quasiparticle operators. In matrix notation it is given by 

with F being a unitary 2Dx2D dimensional matrix. The vacuum \F) for a set of 
quasiparticles obtained by the HFB transformation F is defined in the usual way 
by the request that application of any quasiparticle annihilation operator onto it 
should yield zero, 

a„(F)|F) = Va = !,...,£>. (3) 
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An 77,-quasiparticle state with respect to the vacuum is defined by 

\F{a%)=(flai{F)]\F) for n = l,...,D. (4) 



v«=i 



The vacuum as well as the n-quasiparticle states contain components of many nu- 
cleon configurations with different angular momenta, different angular momentum 



z-components, both parities and various proton and neutron numbers |jT^. In order 
to get physical states which are characterized by good quantum numbers, projection 
techniques have to be used for selecting the components with the desired symmetries. 
The corresponding projection operator Qj^^ is a product of projection operators, 

QiJf ^ P{IM- i^)P(2T,)P(A)P(vr), (5) 

where the first projects on good angular momentum I with 2;-component M, the 
next ones on good isospin ^-component T^, good mass number A and finally on good 



parity tt. Explicit expressions for them can be found in [Q. All projectors apart 
from P(7r) are integral operators. 

Physical configurations are obtained by applying the operator QfJ^^^ on the above 
introduced quasiparticle configurations. Since any of the n-quasiparticle configura- 
tions |F{a^}„) built on the vacuum \F) can be looked upon at the same time as 
a vacuum for another HFB transformation F', it is enough to consider only HFB 
vacua in the following. Thus, a physical configuration with a good symmetry S, 
where S represents the quantum numbers AT^I'", is given by 

\F-SM)= Y: ^MK\F)fl (6) 

K=~I 

The sum over all intrinsic angular momentum ^-components K has to be taken in 
order to avoid a dependency of the projected wave function on the orientation of 
the intrinsic reference frame 0. 

In the VAMPIR approach one symmetry-projected vacuum (^) is used to describe 
the energetically lowest state of a given symmetry S. The configuration mixing 
degrees of freedom and the underlying HFB transformation F are determined 
by variation 

s. AF-^SM\H\F,SM) 
=^ {F;SM\F;SM) = 

This variation leads to three sets of equations which have to be solved self-consistent- 
ly [0 . One obtains the optimal description of the considered yrast state which can 



be achieved by using only a single determinant. Generally this VAMPIR approach 
yields already a rather good approximation to the actual state 
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A similar approach, the FED VAMPIR, uses a hnear combination of a few symmetry- 
projected vacua to describe one yrast state. Here, the correlating configurations, 
which contribute to one state, are searched for by successive variational calculations 
O. Furthermore, by introducing orthogonality constraints, the VAMPIR and the 
FED VAMPIR procedures can be easily extended to the description of excited states. 
These EXCITED VAMPIR and EXCITED FED VAMPIR models can 



describe states of any arbitrary structure. But since each state is searched for by 
variation, a numerical application is quite time consuming. Thus the VAMPIR 
models are especially well suited for the description of the lowest few states of a 
given symmetry S. 



For problems which require a complete excitation spectrum with respect to some 
transition operator, it is better to follow another avenue. Since here only specific 
configurations are needed, one can describe them by expanding the nuclear wave 
functions around a suitable VAMPIR vacuum. This is explained in the next section. 



2.2 The MONSTER on VAMPIR approach 



For states with symmetry S the configuration space is chosen to consist of the 
symmetry-projected vacuum 

\F;SMK)^Qij,\F) (8) 

and the symmetry-projected two-quasiparticle states with respect to it 

\Fa/3; SMK) ^ Ql,^ai{F)a\{F)\F) . (9) 

The underlying HFB transformation F is fixed in a preceding VAMPIR calculation 
for an yrast state, which not necessarily needs to have just the same symmetry S. 
In many cases it is sufficient to use the HFB transformation obtained for the ground 
state of the considered nucleus. 



A general wavefunction for excited states is then given by 
mF)-SM) = \ \F;SMK)g^^.^,+ ^ \FaP; SMK)g'^^A . (10) 

K=-I,...,+I a</3 
^, K=-I,...,+I J 

The expansion coefficients are obtained by diagonalizing the effective Hamiltonian 
in the space of the nonorthogonal configurations (||, |^). 

The above choice of the model space restricts the MONSTER on VAMPIR approach 
to excited states with a structure similar to that of the underlying projected vacuum. 
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This is here desired, since we are interested in states predominantly populated by 
a transition from the vacuum via a one-body operator. Furthermore it should be 
mentioned that the MONSTER approach provides an approximate elimination of 
spurious center-of-mass excitations, which can be introduced by the configuration 
mixing if more than one major oscillator shell is taken as particle basis (for details 
see 0). 



2.3 Symmetry restrictions imposed on the 
HFB transformation 



If the most general HFB transformations are allowed, symmetry-projected vacua of 
type (|^) can describe arbitrary states in arbitrary nuclei |T^. Up to now, however. 



this has not been achieved in any numerical implementation. Instead, for the existing 
numerical realizations of the VAMPIR approaches, certain symmetry requirements 
were imposed on the underlying HFB transformations. 



In the first VAMPIR calculations ||T0[ axial symmetry and time-reversal invariance 
were assumed, parity and proton-neutron mixing were neglected and only real HFB 
transformation coefficients were admitted. As a consequence of these approximations 
the real VAMPIR approach was only suitable for positive parity states with even 
angular momenta in doubly even nuclei. Performing MONSTER calculations on 
top of such a real VAMPIR solution these limitations are removed. In the two- 
quasiparticle approximation not only states with arbitrary spin-parity in doubly 
even but also in doubly odd systems become accessible. However, for the description 
of odd spin and/or negative parity states in doubly even nuclei one has to use mean 
fields obtained for different spin and/or parity values than the considered one. And 
for doubly odd systems one even has to rely on mean fields for neighbouring systems. 
Considering as configuration space for the MONSTER calculation all symmetry- 
projected one-quasiparticle states 

\Fa;SMK)^Qf,^aUF)\F), (11) 

arbitrary states in odd systems become approachable, too. Here again, one has to 
build the one-quasiparticle spectrum on the mean field determined for a neighbouring 
doubly even nucleus. 

In the more recent implementations of the VAMPIR calculations parity as well as 
proton-neutron mixing are taken into account and essentially complex HFB trans- 
formations are admitted. In this complex VAMPIR approach only axial symmetry 
and time reversal invariance are kept |T^. This introduces many more correlations 



into the projected vacua than the older real calculations. Now states of arbitrary 
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spin-parity in both doubly even and doubly odd systems are accessible. The com- 
plex vacua contain all possible two nucleon couplings, but some four-, six- and more 
nucleon couplings are missing [0. Only recently the MONSTER approach could 
be extended to such complex VAMPIR vacua [|T^ . Here, in addition to doubly even 
nuclei, also the calculation for doubly odd nuclei can be based on a HFB transforma- 
tion determined for just the considered system itself. And the structures which are 
missing in the complex VAMPIR vacuum can be introduced by the two-quasiparticle 
admixtures. However, for the description of odd systems one still has to use the one- 
quasiparticle approximation like in the older real approach. But now aside from the 
mean fields of neighbouring even-even nuclei also the ones of neighbouring odd-odd 
nuclei can be chosen to construct the one-quasiparticle states (0). 

Generally, since the underlying HFB transformation is much more general and the 
configuration space is much larger, the complex MONSTER is expected to yield a 
better description of nuclear spectra than the real one. 



3 A systematic application to sd shell nuclei 



In [0 we presented as a first test of the new complex MONSTER approach an 
application to ^°Ne and ^^Ne. Considering the sd shell as model space, it was 
shown that the model is able to reproduce the shell model spectrum of ^°Ne exactly. 
For ^^Ne it yielded an excellent approximation to the shell model result. A clear 
improvement on the older real approach was found. However, the two investigated 
nuclei are even-even systems. Moreover, they are quite light : for both nuclei the 
number of shell model states with a certain spin is comparable to the number of 
complex MONSTER configurations. How does the approach perform for heavier 
nuclei, where the number of states is much larger than the number of complex 
MONSTER configurations ? Furthermore, which kind of agreement is obtained 
for odd-odd nuclei and for odd mass systems, where the simpler one-quasiparticle 
description is used ? These questions are adressed in the present investigation. 

Again the IsOd shell was chosen as single particle basis in order to be able to compare 
with exact shell model configuration mixing calculations. The single particle energies 
e {d5/2) = -4.15 MeV, e (sl/2) = -3.28 MeV, and e {d3/2) = +0.93 MeV for both 
protons and neutrons |T8[ and the effective two-body residual interaction (the mass- 
dependent version of the Chung and Wildenthal force ||T9| {V{A) = V{18) x (18/^4)") 
with a = 1/3 instead of a = 0.3) are the same as in our earlier calculations. 
This force is generally accepted to be "the standard" force for the sd model space. 
The results are compared in addition to those obtained by the more restricted real 
MONSTER approach. 
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3.1 Even-even nuclei 



As examples for "heavier" even-even sd-shell nuclei, we calculated the spectra of 
l^Mg, llMg, llSi and ^^Si within the complex and the real MONSTER approaches. 
For each MONSTER calculation the underlying HFB transformation was deter- 
mined by a preceding VAMPIR calculation for the 0+-ground state. This HFB 
transformation F(0+) has then been used for states of all spins in one nucleus. 



The energies obtained by the complex and the more restricted real MONSTER have 
been compared to those from an exact shell model diagonalization. The result of 
this investigation is summarized in table 1. The average deviation of the energies 
of the yrast states with spin 0+ to 5+ from the shell model energies is presented 
for each nucleus. The deviation is given in keV and in percent of the shell model 
ground state energy of the nucleus referred to. The exact ground state energy is 
also displayed. For comparison the corresponding values for the nuclei ^''Ne and 
^^Ne of our first application are listed here, too. The table shows furthermore the 
minimum and maximum number of shell model configurations per spin for each 
nucleus and the corresponding spin value. In addition the number of shell model 
states per spin averaged over the states with spin 0+ to 5"*" is given. It should be 
compared to the number of MONSTER configurations available in the S(i-shell. As 
explained in more detail in |TB[ the complex MONSTER provides 57 configurations 
for the description of spin O"*" states, 151 for 1+, 223 for 2+, 259 for 3+, 275 for 4+, 
and always 277 for all higher spin states. In the more restricted real MONSTER 
approach less configurations per spin are available : 21 for the / = states, 30 for 
/ = 1, 61 for / = 2, 56 for / = 3, always 60 for / = 5, 7, 9, 11, 13, and always 73 for 
the / = 4, 6, 8, 10, 12, 14 states, respectively. 



As expected for both MONSTER approaches the mean deviation of the yrast state 
energies from the exact shell model result, averaged for spins O"*" to 5"*", does increase 
with the number of actually existing shell model states. For the complex MONSTER 
the largest deviation measured relatively to the shell model ground state energy 
occurs for ^^Mg and amounts to 1.43% . For ^°Ne the complex calculation reproduced 
the exact result. In the real case the largest average deviation, 2.42 % , is found 
for 24Mg and the smallest, 0.87% , for '^'^Ne. For all nuclei the complex MONSTER 
approach yields a clear improvement on the real one. The largest energy gain could 
be achieved in ^^Mg, where the averaged deviation of the real MONSTER is 3.8 
times larger than the one of complex MONSTER. The smallest amount of energy 
was won for ^^Mg, where the average deviation of the real approach is by a factor 
of 1.3 larger compared to the complex one. For all even-even nuclei the complex 
MONSTER approach yields an excellent approximation to the exact calculation. 

As an example, in fig. 1, the spectra of ^^Si as obtained by the complex and the real 
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MONSTER are shown up to ~ 14 MeV above the shell model ground state energy. 
For comparison the energies resulting from the shell model diagonalization are plot- 
ted up to ~ 12 MeV excitation energy. It can be seen that the complex MONSTER 
reproduces the shell model spectrum well and shows a clear improvement on the real 
MONSTER : the states are considerably more bound and many more shell model 
states can be described. 



The nucleus ^^Si was investigated in more detail, since it has the largest number of 
states of all S(i-shell nuclei. There exist in total 93710 states. More precisely, there 
are 3372 r=0+ states, 9216 r=l+, 13562 /^=2+, 15385 r=3+, 15089 r=4+, 
12876 1'"= 5"^ and 9900 J'^= 6^ states. For the higher spins the number of states 
per spin continues to decrease; there is just one state with the highest possible spin 
I'^ = 14"*". In fig. 2 the yrast state energies as obtained by various models of the 
VAMPIR and MONSTER family are compared for the spins 0+ to 6"^. The energies 
from the exact shell model calculation are given in the leftmost column. Then, 
from left to right, the energies calculated by the complex MONSTER, the complex 
VAMPIR, the real MONSTER, and the real VAMPIR are shown. Note that the 
real VAMPIR approach can only describe even spin states. 

The average deviation of the energies obtained by the complex MONSTER approach 
from the shell model, now including the 6"'" yrast state, is about 1085 keV or 0.79% 
of the shell model ground state energy of 137.495 MeV. The deviation for the even 
spin states is on average somewhat smaller, namely 840 keV, than the one for the 
odd spin states, which amounts to 1412 keV. This was found also for the other 
investigated even-even nuclei, apart from ^^Mg, where the difference is insignificant. 
There the average deviation of the odd spin states was only 51 keV larger than 
for the even spin states. Real MONSTER gives for ^®Si an average deviation of 
1746 keV, which amounts to 1.27% of the shell model ground state energy. 



As it is to be expected, the complex MONSTER yrast states are more bound than 
the ones calculated by complex VAMPIR, except for the 0"*" VAMPIR state, which 
is stable with respect to any projected two-quasiparticle excitations of the same 
symmetry. Obviously the odd spin states l"*", 3"*" and 5"^ are not well described by the 
complex VAMPIR. When the difference between the binding energies obtained by 
the two models is so pronounced, it can be inferred that these state are dominated 
by the above mentioned "missing couplings", which are not accessible within the 
complex VAMPIR model. This was also found in Ref. [16|, where the odd spin yrast 
states of ^^Ne were rather poorly described by the complex VAMPIR approach. 



The differences in the binding energies of real MONSTER and real VAMPIR are 
negligibly small. The even spin real MONSTER yrast states, 2+, 4+ and 6^, are 
on average only 56 keV more bound. The O"*" state does not mix with the projected 
two-quasiparticle states of the same symmetry in the real case, too. 
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3.2 Odd-odd nuclei 



In order to study the quality of the complex MONSTER approach for apphcations 
to odd-odd nuclei, we calculated the energy spectra of ^gF, ^^Na, If'Na, 13AI, ^gAl 
and f^P. For each nucleus the complex MONSTER configuration space was built on 
the HFB transformation, which yielded the energetically lowest VAMPIR solution 
in just this nucleus, no matter which spin it corresponded to. The spectra have 
also been calculated with real MONSTER. In this case one has to take the HFB 
transformation of a neighbouring even-even nucleus to construct the MONSTER 
configuration space. Naturally, one has several possibilities : for an odd-odd nucleus 
with Z protons and N neutrons one can take the even-even neighbours with the 
proton and neutron numbers (Z-1,N-1), (Z-1,N+1), (Z+1,N-1), and (Z+1,N+1). 
It is a priori not clear, which choice may be the best. Therefore we looked in 
each neighbour for the HFB transformation, which belongs to the 0"*" ground state 
VAMPIR solution and used each for a subsequent real MONSTER calculation. 

The energies obtained by both MONSTER approaches have been compared to the 
ones of a full shell model configuration mixing calculation. To get a quantitative 
measure, for each nucleus as already for the even-even nuclei, the average deviation of 
the energies of the yrast states with spins O"*" to S"*" from the shell model energies has 
been calculated. The resulting deviations are shown in table 2, again in keV and in 
percent of the shell model ground state energy. Besides the energy the table displays 
the spin of the shell model ground state and, for complex MONSTER the spin of 
the VAMPIR solution which was used. For real MONSTER always the result for 
that spectrum which yielded the best agreement with the shell model is listed in the 
table. The corresponding neighbouring nucleus is indicated. The real MONSTER 
spectra for one odd-odd nucleus, obtained from the HFB transformations of the 
various possible neighbours, differ often quite much. So, e.g., for ^^Na the average 
deviation of the yrast states 0^ to 5^ is 4.85% if the HFB transformation for the O"*" 
VAMPIR solution of ^^Ne is used, and 1.47% for the one of ^^Ne. 

Table 2 also shows the minimum and maximum number of shell model configurations 
per spin and the number per spin averaged for spins 0+ to Since protons and 
neutrons are mixed by the complex HFB transformation, the same configurations are 
used for the calculation of even-even and odd-odd nuclei in a complex MONSTER 
calculation. Thus the number of configurations contributing to one spin given in 
section ( |3.1| ) is valid here, too. For the real MONSTER with no proton-neutron 
mixing the situation is different compared to the even-even case. For the description 
of odd-odd nuclei in sd shell there are 14 configurations available for the 0+ states, 
39 for 1+, 57 for 2+, 67 for 3^, 71 for 4+ and always 72 for the states with spins 
5"*", . . . 13"''. These numbers are comparable in size to the ones of the even-even case. 
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As already for the even-even nuclei, also for the odd-odd ones an increase of the 
average deviation with the number of shell model states per spin is observed. For 
the whole shell model spectrum is exactly reproduced by complex MONSTER 
using any complex HFB transformation F. This is to be expected, since ^°F belongs 
to the same isospin multiplet like ^''Ne. ^°Ne encloses states of all isospins T = 0, 1, 2, 
^°F only the ones with isospins T = 1 and 2. This provides another stringent test 
for the numerics. The largest deviation relative to the shell model ground state 
energy was found for ^^Na : 2.21%. ^^Na is a typical example for which one should 
better choose different HFB transformations for states with different spins, instead 
of using one transformation for all states. If on uses the HFB transformation of the 
O"*" VAMPIR solution for the low spin states 0"^ to S"*" and the transformation F(4"'") 
for the higher spins, the average deviation of the yrast states 0+ to S"*" can be reduced 
to 995 keV or 1.26% of the shell model ground state energy. The possibility to choose 
different HFB transformations for different spins should always be considered, since 
we found that the complex MONSTER energies are quite sensitive on the underlying 
HFB transformation in contrast to the real case [|16|. 

The average deviations of the real MONSTER are always bigger than for the com- 
plex approach, apart from ^^Na, where both are of similar size. But this could be 
improved by using several HFB transformations in the complex case. As far as the 
higher excited states are concerned, the complex MONSTER is always superior to 
the real one, since it is able to describe many more states. 

Compared to the even-even case the deviations obtained by complex MONSTER 
are somewhat bigger for each mass number. It is known that the structure of 
odd-odd nuclei is generally more complicated. More types of couplings contribute 
significantly. This seems to be reflected in our results. But also for the odd-odd case 
the exact spectra are over all rather well approximated by the MONSTER approach. 

In fig. 3 the energy spectra of ^^Na calculated with the complex and the real MON- 
STER are plotted up to spin 7 and an energy of —52 MeV. For comparison the 
lowest shell model energies for each spin are presented. In the real case, MONSTER 
is built on the HFB transformation determined for the O"*" VAMPIR solution in ^''Ne. 
This gave the best agreement with the shell model. The most bound complex VAM- 
PIR solution in this nucleus is the one for the 1"^ yrast state. Complex MONSTER 
on top of the corresponding HFB transformation corrects then the sequence of the 
yrast states to the one of the exact calculation, where the ground state has spin 3~^. 
The largest deviation of a complex MONSTER yrast state energy from the shell 
model value occurs for the O"*" yrast state : 362 keV, the smallest for : 83 keV. 
Complex MONSTER shows also for the higher lying states an excellent agreement 
with the energies of the exact shell model and yields a significant improvement on 
the real MONSTER. 
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Finally, as an example for the heavier odd-odd nuclei, fig. 4 displays the calculated 
energies of ^°P. In this figure two real MONSTER calculations are presented, the 
ones yielding the best and the worst agreement with the result of the exact diago- 
nalization. The former has been got by using the HFB transformation for the O"*" 
VAMPIR solution in ^^Si, the latter by the HFB transformation for the 0+ VAMPIR 
solution in ^^S. The lowest complex VAMPIR solution in ^°P was found to be the 
one with spin O"*". Its HFB transformation was used in the complex calculation. 
The complex and the real MONSTER results are plotted up to ~11 MeV above the 
shell model ground state for spins 0"'' to 5~^. For each spin the lowest energies of 
the exact diagonalization are given. Complex MONSTER reproduces the general 
trends of the spectrum rather well. The largest deviation of the yrast state energies 
is found for spin 2"*" and amounts to 1579 keV or 1.01% of the shell model ground 
state energy, the smallest, 1157 keV or 0.74%, for spin O"*". The complex calculation 
shows a clear improvement on both real ones, where much less states are described. 



3.3 Odd mass nuclei 

As described in section ( |2.3| ), in the case of odd mass nuclei the complex MONSTER 
spectrum has to be built on a HFB transformation determined for a neighbouring 
even mass nucleus. We proceeded here in the following way. We calculated the 
shell model spectra of several of those odd mass nuclei, which are neighbouring 
the even mass nuclei discussed in sections ( |3.1] ) and Then we used the HFB 

transformations, which we had already determined, to built up the MONSTER 
one-quasiparticle configuration spaces. Accordingly, for the mass parameter of the 
Chung- Wildenthal force, we took always the mass of the even mass nucleus of which 
the HFB transformation was used. The energy spectra of ijNe, ifNa, fl^S^ is^l 
and flSi have been calculated. For ^^Ne the HFB transformations of the deepest 
VAMPIR solutions in ^^Ne and ^^Na have been used. The other odd mass nuclei 
with mass number A have been calculated always in four different ways : using the 
HFB transformation for the deepest VAMPIR solution of the even-even neighbour 
with mass A— 1, the one of the even-even neighbour with mass A+1 and the ones 
of the odd-odd neighbours with A— 1 and A+1. Thus, e.g., ^^Si has been calculated 
using the mean fields of ^^Si, ^^Si, ^'^Al and ^^P. In the case of real MONSTER only 
the mean fields of the two neighbouring even-even nuclei are available. Both have 
been used for the real calculations. 

The MONSTER configuration space for the calculation of odd mass nuclei is much 
smaller than the one for the even mass nuclei, because there are much less one- 
quasiparticle than two-quasiparticle excitations. For complex MONSTER in the sd 
shell there are 12 configurations available for states with spin 1 = 1/2, 20 for 1 = 3/2 
and always 24 for states with higher spins. The real MONSTER provides even only 
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3 configurations for 1—1/2 states, 5 for 7 = 3/2, and always 6 for all higher spin 
states. 

Let us first turn our attention to ^^Ne. It is the lightest odd mass nucleus which we 
have calculated. For this nucleus there are in total 1935 shell model configurations. 
More precisely, there are 199 1= 1/2, 341 / = 3/2, 400 / = 5/2, 368 / = 7/2, 287 
1 = 9/2, 183 1 = 11/2, 100 1 = 13/2, 41 7 = 15/2, 14 7 = 17/2, and 2 7 = 19/2 states. 
Fig. 5 shows for each spin the lowest energies resulting from the complete shell 
model diagonalization up to spin 7= 11/2. In addition the results of two complex 
MONSTER calculations are displayed up to an energy of —38 MeV. One calculation 
was done by using the mean field of the 0"*" ground state for ^^Ne, the other is built 
on the mean field of the deepest VAMPIR solution for ^^Na, i.e., the 1+ solution. 
In fig. 5 also the result of a real MONSTER calculation is plotted for comparison. 
It is built on the HFB transformation determined for the 0+ ground state of ^^Ne 
by a reai VAMPIR calculation. The obtained agreement between both complex 
MONSTER results and the shell model is very good. For the calculation built on 
^^Ne the average deviation of the yrast state energies from the exact ones is 490 keV. 
That is 1.03% of the shell model ground state energy, which is —47.753 MeV. The 
largest deviation occurs for the 11/2 state and amounts to 664 keV or 1.39% of the 
shell model ground state energy, the smallest for spin 5/2, and is 256 keV or 0.54%. 
The calculation using the mean field of ^^Na gives an average deviation of 463 keV 
or 0.97% measured relative to the shell model ground state energy. The largest 
deviation, 653 keV or 1.37%, was found for the 1/2 state, the smallest, 308 keV or 
0.64%, for the 7/2 yrast state. Thus, the mean field of the 1+ yrast state of ^^Na is on 
average shghtly better suited for the description of the yrast states of ^■'^Ne than that 
of the 0"*" ground state of ^^Ne. Except for the spin 1/2 and 5/2 yrast states, it yields 
deeper bound states. Also the excited states obtained by this calculation display 
the bigger binding energies. The real MONSTER gives for the yrast state energies 
an average deviation of 576 keV or 1.21% from the shell model ones. The 1/2 state 
shows the largest deviation : 704 keV or 1.47%, the 9/2 state the smallest : 467 keV 
or 0.98%. For most of the yrast states both complex MONSTER calculations yield 
more binding than the real ones, as it is expected. However, e.g., for the 9/2 yrast 
state the real MONSTER yielded the deepest energy. The complex results, built on 
the mean fields of ^^Na and ^^Ne, are higher by 16 keV and 114 keV, respectively. 
This indicates that in the complex case, the mean field determined for a low spin 
may not be well suited for the description of high spin states. 

As example for an heavier odd mass nucleus, we present the energy spectra of ^^Si 
obtained by using the mean fields of the mass 28 systems. The number of shell model 
states in this nucleus is huge compared to the number of MONSTER configurations 
available. In total this nucleus has 80112 shell model states. 5638 of these states 
have spin 7 = 1/2, 10176 7 = 3/2, 12877 7 = 5/2, 13450 7 = 7/2, 12240 7 = 9/2, 
9835 7 = 11/2, 7053 7 = 13/2, 4469 7 = 15/2, 2502 7 = 17/2, 1197 7 = 19/2, 485 
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I = 2112, 152 7 = 23/2, 35 7 = 25/2, and 3 are spin I = 27 12 states. For each spin 
up to spin 9/2 the lowest energies resulting from the exact shell model calculation 
are shown in fig. 6. The obtained MONSTER spectra are plotted up to an energy 
of ~ 10 MeV above the shell model ground state energy of —147.020 MeV. One 
complex MONSTER spectrum is obtained by using the HFB transformation of the 
VAMPIR solution for the O"*" ground state in 28Si, the other uses the transformation 
of the O"*" VAMPIR solution in ^^Al. The calculation on top of the mean field of the 
even-even nucleus shows an average deviation of the yrast state energies from the 
exact ones of 1373 keV or 0.93% in terms of percentage of the shell model ground 
state energy. The yrast state with spin 9/2 shows the largest deviation, 1575 keV or 
1.07%, the 7/2 yrast state the smallest, 1291 keV or 0.88%. The average deviation 
obtained by using the mean field of the odd-odd nucleus is somewhat bigger, namely 
1912 keV or 1.30%. The largest and smallest deviations, which occur, are 2172 keV 
or 1.48% and 1607 keV or 1.09%, respectively. The general trend of the shell model 
spectrum is reproduced by both complex calculations. In fig. 6 in addition also a real 
MONSTER spectrum is shown. It was obtained by using the HFB transformation 
of the O"*" VAMPIR solution in ^*Si. This spectrum yielded a better agreement with 
the shell model result than the real calculation on top of the mean field for ^°Si. 
For the displayed real MONSTER spectrum the average deviation of the yrast state 
energies from the exact ones is 1973 keV or 1.34%. The largest deviation occurs for 
the 1/2 yrast state, 2286 keV or 1.55%, the smallest for 7/2, 1640 keV or 1.12%. 
The complex MONSTER calculation on top of the mean field of ^^Si shows a clear 
improvement over this best real result. 



For all the considered odd mass nuclei, apart from ^^Ne, we found that by using 
the HFB transformation of a neighbouring even-even nucleus, a better agreement 
between the complex MONSTER result and the shell model could be achieved than 
by using the mean field of an odd-odd neighbour. In general the odd mass spectra 
on top of a mean field of a neighbouring even-even nucleus approximated the shell 
model results well (more details can be found in 0]). The only exceptions were the 
spectrum of ^^Na obtained by using the HFB transformation of the 0"*" VAMPIR 
solution for ^^Ne and the spectrum of ^^Mg with the transformation for ^^Mg. These 
two calculations gave a rather poor agreement with the exact result. In both cases 
it could be much improved by using the HFB transformation of ^^Mg. Thus, in the 
odd mass case the choice of the neighbour can be crucial. One should always try 
several ones to find the mean field, of which the structure is best suited for the odd 
mass nucleus under consideration. 



Finally we investigated for the odd mass case the dependence of the complex MON- 
STER energies on the spin of the mean field determined by the preceding VAMPIR 
calculation. For ^^Ne we built always the full MONSTER spectrum on each of the 
yrast solutions of ^^Ne obtained by complex VAMPIR. The same was done for real 
MONSTER. Naturally here only mean fields for even spin states are available. 
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Fig. 7 shows a typical result of this investigation. The energies of the five lowest 
3/2 states of ^^Ne are plotted in dependence on the spin used in the preceding 
VAMPIR calculation for ^SNe. For comparison, the leftmost column presents the 
energies of the three lowest 3/2 shell model states. Solid lines denote the results of 
the complex MONSTER, dotted ones those of the real counterpart. It can be seen 
that real MONSTER displays only a weak dependence on the spin of the underlying 
HFB transformations. For the complex MONSTER this is also true if the VAMPIR 
transformations of the low spin mean fields, O"*", 2"*" and 4"*", are considered. The 
mean fields of the higher spin states are insufficient for the description of the 3/2 
states. The poor description of the 3/2 states using the VAMPIR transformations 
of the low odd spins, 1"*", 3"'' and 5"*", is due to the fact that those yrast states are 
dominated by structures which are missing in a complex VAMPIR calculation . 

The same behaviour as for the 3/2 states was also found for most of the spins of 
^^Ne. Only for the very high spins, where the number of MONSTER configuration 
is of similar size as the number of shell model states, the dependence becomes small. 
When the number of MONSTER configurations is bigger than the number of shell 
model states, the resulting energies are independent of the underlying transforma- 
tion. 

This investigation confirms what we have found already for the even mass case . 
The complex HFB transformation shows a much stronger dependence on the spin of 
the underlying VAMPIR calculation than the real one. It means that the structure 
of the complex mean field differs more from spin to spin. The choice of the spin of 
the underlying VAMPIR calculation is important in the complex description of odd 
mass nuclei, too. 



4 Summary 

A recently developed nuclear structure approach, the complex MONSTER, has been 
studied. The approach is the newest member of a group of models, which are 
based on variational methods. This group can be divided into two subgroups, one 
consisting of the VAMPIR models, the other of the MONSTER approaches. Both 
are designed for large scale nuclear structure calculations. 

The models of the VAMPIR family use symmetry-projected HFB quasiparticle vacua 
as test wavef unctions. The underlying HFB transformations are determined by vari- 
ation. By construction, these models can only be used to describe the lowest few 
states of a certain symmetry. In the first numerically realized VAMPIR models 
the HFB transformations were rather restricted : time-reversal invariance and ax- 
ial symmetry were required, parity and proton-neutron mixing were neglected and 
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only real HFB transformation coefficients were allowed. In the newest models most 
restrictions are removed and essentially complex HFB transformations are allowed. 
Only time reversal and axiality are kept. 

In the MONSTER approaches the nuclear wavefunctions are expanded around a 
VAMPIR solution for the ground or an yrast state. The spectrum of excited states 
is obtained by diagonalizing the Hamiltonian in the space of the VAMPIR solu- 
tion and all symmetry-projected two-quasiparticle configurations with respect to it. 
These models are therefore suited for problems, where a complete set of excitations 
with respect to a particular transition operator is needed. Till lately the MON- 
STER approach was restricted to the use of real VAMPIR solutions. In a recent 
work it was generalized. Now complex HFB transformations can be used to build 
up the configuration space. Then more nucleon correlations are considered already 
in the mean field and the configuration space becomes much larger. Consequently 
a more detailed description of nuclear spectra is expected by this complex MON- 
STER approach. Furthermore one can use HFB transformations determined for 
the particular system under consideration not only for even-even systems, as in the 
real approach, but also for odd-odd systems. Only for odd mass nuclei, which are 
described in the one-quasiparticle approximation, one still has to use the mean field 
of a neighbouring nucleus, which now may be even-even or odd-odd. To get insight 
into the ability of this complex approach a systematic investigation was performed 
in this work. 

We presented first results of an application of the complex MONSTER approach 
to several nuclei in sd shell. This basis, which is rather tiny for a MONSTER 
calculation, has been chosen to enable a comparison with complete shell model 
configuration mixing calculations. The investigated even-even and odd-odd nuclei 
have masses between A=20 and 30, the odd mass nuclei reach from A=21 to 29. 
Thus, quite light systems with only a few valence nucleons are considered as well as 
heavier systems, up to ^^Si with the biggest number of shell model configurations in 
sd shell. 

For the even-even nuclei the complex MONSTER yielded an excellent agreement 
with the exact shell model approach. The largest mean deviation of the yrast state 
energies from the shell model result, averaged over the states with spins to 5, was 
found for ^^Mg and amounts only to 1.43% of its shell model ground states energy. 
Concerning the odd-odd systems it was demonstrated that complex MONSTER 
reproduced the shell model spectrum of ^°F exactly. ^^Na showed the largest average 
deviation of the energies of the yrast states with spins to 5 from the shell model 
ones, namely 2.21% of the ground state energy. The average deviations of the odd- 
odd systems were slightly bigger than the ones of the even-even systems with the 
same mass number. But in general the complex MONSTER yielded a very good 
approximation to the shell model result also for odd-odd nuclei. For both even-even 
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and odd-odd nuclei a clear improvement on the previous real approach was proved. 

For odd mass nuclei it was found that the description depends much on the un- 
derlying mean field. By choosing a suitable mean field a good agreement with the 
shell model could be achieved for all investigated nuclei. As examples the results of 
^^Ne and ^^Si were presented. For ^^Ne, using the HFB transformation of the l"*" 
VAMPIR solution in ^^Na, the deviation of the yrast states energies, averaged over 
spins 1/2 to 11/2, was 0.97% of its shell model ground state energy. For ^^Si, with 
the mean field determined for the ground state of ^^Si, the average deviation of the 
yrast states energies for spins 1/2 to 9/2 amounted to 0.93% relative to the ground 
state energy. Again the complex approach achieved better agreement with the shell 
model than the real one. 



Already the previous real MONSTER has been applied successfully to the descrip- 
tion of nuclear structure phenomena in various mass regions, e.g., in the mass 130 
region 
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The present investigation has clearly shown that the new complex 
approach is even superior. So we expect that the complex MONSTER approach will 
develop as a more powerful tool for nuclear structure studies in large model spaces. 



Up to now only the energies calculated by the different models have been compared. 
It will be interesting to look also for other observables, like transition strengths up to 
high excitation energies, and to analyze the wavefunctions. This is planned for the 
future. Finally we would like to mention that the description of nuclear wavefunc- 
tions by VAMPIR and MONSTER could be improved by removing all restrictions 
imposed on the HFB transformation. The development of such approaches allowing 
the most general HFB transformations is in progress, too. 



We thank Prof. Dr. Herbert Miither for performing the shell model calculations with 
the Glasgow code. Furthermore this work was partly supported by the Graduier- 
tenkolleg "Struktur und Wechselwirkung von Hadronen und Kernen" (DFG, Mu 
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Figure Captions 



Fig. 1 Comparison of the energy spectra of ^^Si for spins to 6 as obtained by three 
different approaches : the shell model configuration mixing approach (SCM), 
the complex MONSTER (CM) and the real MONSTER (RM). The energy is 
given relative to the ^^O core. For each spin only the lowest few shell model 
states arc presented. The MONSTER results are displayed up to excitation 
energies of ~ 14 MeV above the shell model ground state. The complex and 
the real MONSTER calculation has been built on the corresponding complex 
and real VAMPIR solution for the O"*" ground state. 

Fig. 2 The energies of the yrast states with spin 0^ to 6^ in ^^Si as obtained by 
various approaches : the shell model (SCM), the complex MONSTER (CM) 
(based on the 0+ solution of complex YAMPIR), the complex YAMPIR (CV), 
the real MONSTER (RM) (on top of the 0+ solution of real VAMPIR), and 
the real VAMPIR (RV). In the last approach only the even spin states are 
accessible. Spin and parity are indicated on the l.h.s. of each level. The 
energy is given relative to the ^^O core. 

Fig. 3 Same as Fig. 1, but for the odd-odd nucleus ^^Na and up to spin 7. Each 
MONSTER spectrum is shown up to 7.5 MeV above the shell model ground 
state. The complex MONSTER was built on the complex VAMPIR solution 
obtained for the 1+ yrast state, the real MONSTER was based on the real 
VAMPIR solution for the 0+ ground state of ^^Ne. 

Fig. 4 Same as Fig. 1, but for the odd-odd nucleus ^°P and up to spin 5. The 
results of two real MONSTER calculations are displayed. Each MONSTER 
spectrum is shown up to 11 MeV above the shell model ground state. The 
complex MONSTER was built on the complex VAMPIR solution obtained for 
the 0+ yrast state. One real MONSTER calculation has been based on the 
VAMPIR transformation for the 0+ ground state of ^^Si, the other on the one 
of the 0+ ground state of ^^S. 

Fig. 5 Same as Fig. 1, but for the odd mass nucleus ^^Ne and for spins 1/2 to 
11/2. The results of two complex MONSTER calculations are displayed. Each 
MONSTER spectrum is shown up to -38 MeV. One complex MONSTER 
calculation has been based on the complex VAMPIR transformation for the 
O"*" ground state of ^^Ne, the other was using the one for the l"*" yrast state of 
^^Na. The real MONSTER calculation was built on the real VAMPIR solution 
for the 0+ ground state of ^^Ne. 

Fig. 6 Same as Fig. 1, but for the odd mass nucleus ^^Si and for spins 1/2 to 
9/2. The results of two complex MONSTER calculations are displayed. Each 
MONSTER spectrum is shown up to 10 MeV above the shell model ground 
state. One complex MONSTER calculation has been based on the complex 
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VAMPIR transformation for the O"*" ground state of ^^Si, the other was using 
the one for the O"*" yrast state of ^^Al. The real MONSTER calculation was 
built on the real VAMPIR solution for the 0^ ground state of ^^Si. 

Fig. 7 The five lowest 3/2+ MONSTER states of ^^Ne are plotted versus the spin 
of the yrast state in ^^Ne, of which the VAMPIR transformation was used in 
each calculation. Solid lines refer to the complex MONSTER, dotted lines to 
the real MONSTER results. In the latter case only transformations for even 
spin values are available. For comparison the lowest three 3/2+ shell model 
(SCM) energies are also given. The energy is given relative to the ^^O core. 
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Table 1: 



A comparison of the deviations of the yrast state energies calculated with complex 
and real MONSTER from the exact shell model energies is provided for the considered 
even-even nuclei. In addition the minimum, the maximum and an average number of 
shell model states are shown. 
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Al: average deviation of the energies of the yrast states O"*" to 5"'' calculated with 
complex (CM) and real (RM) MONSTER from the shell model energies in keV. 

SCMp^: shell model ground state energy in MeV. 

A2: same as Al but now in percentage of the shell model ground state energy. 
Nmin'- minimum number of shell model states states per spin and corresponding spin. 
Nmax'- maximum number of shell model states states per spin and corresponding spin. 



leaver- uumbcr of shcU model states per spin averaged over spins O'*' to 5"''. 
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Table 2: 

Same as table 1 but for the considered odd-odd nuclei. 
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In addition to the shorthand notations of table 1, we have used: 

F: the underlying HFB transformation F is indicated. For the complex MONSTER 
(CM) the spin of the corresponding VAMPIR solution is given. In the real case 
(RM) the neighbouring even-even nucleus is shown, which has been used to de- 
termine the HFB transformation. 

SCMg^: contains in addition to the energy, as in table 1, also the shell model ground 
state spin. 



-51 



-53- 



"5 -55 



CD 

c -57 

LJJ F — _ 22m 

Na 



-59- 



_ ___ SCM:- 

— CM(F(r)):- 

RM (FrNe,00): " 

-^^ \ \ \ \ \ \ \ r 

0' V 2" 3" 4' 5' 6' 7" 

Total angular momentum and parity 



-146 
-148 



-150- 



> 

CD 



D) 
CD 

-154- 



-156- 



-158 



30 p 

SCM 

CM (F(00) 
RM (FrSi,00) 
RM (FC'S,00) 



T 

0^ 



T 



T 



T 



T 



T 



r 2' 3' a; 5' 

Total angular momentum and parity 



-37 



-39- 



-41 - 



>.-43|- 

D) 
Q) 

LU 

-45h - - ^^Ne 

SCM 

-471- CM (F(''Ne,00) 

CM(Ff Na,r)) 
RM (FrNe,00) 

-49^^ \ \ \ \ \ r 

1/2 3/2 5/2 7/2 9/2 11/2 

Total angular momentum 



-136 
-138 



-140- 



>,-142^ 

D) 



-144- 



-146- 



-148 



SCM 

CM (FrSi,00) 
CM (FrAI,00) 
RM (FrSi,00) 



T 



T 



T 



9/2 



1/2 



3/2 5/2 7/2 

Total angular momentum 



3/2' states of '^Ne 



-29 r 
-31 '- 
-33 r 
-35E- 
CD -37P 



>> -39F 
D) 

CD 

-41^ 

LU 



-43 r 

-45 r 
-47 r 



-49 



— I r 

SCM 0' 



r 



T" 

2' 



T 

3' 



T" 

4' 



5' 



6' 



-r 

7' 



8' 



T" 

9' 



-r 

10' 



